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Abstract 

Based on a previously found general class of quantum improved exact solutions 
composed of non-interacting (dust) particles, we model the gravitational collapse of 
stars. As the modeled star collapses a closed apparent 3-horizon is generated due to 
the consideration of quantum effects. The effect of the subsequent emission of Hawking 
radiation related to this horizon is taken into consideration. Our computations lead us 
to argue that a total evaporation could be reached. The inferred global picture of the 
spacetime corresponding to gravitational collapse is devoid of both event horizons and 
shell-focusing singularities. As a consequence, there is no information paradox and no 
need of firewalls. 
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1. Introduction 

The evolution of a collapsing star has been one of the favorite research topics in theoret¬ 
ical physics for almost a century. The hrst works in the framework of General Relativity 
(see, for instance, imziisi) seemed to show that, under the appropriate circumstances, the 
formation of a black hole was unavoidable. Moreover, with the later development of the 
singularity theorems [1] it was also conjectured that singularities could exist hidden inside 
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black holes. The contribution of Quantum Theory to this picture was twofold. On the one 
hand, Hawking showed that the event horizon of a black hole instead of behaving as a one 
way membrane could emit radiation [Sj. This opened the possibility that black holes could 
evaporate completely. However, this also created a new problem: What happens to all the 
information that had been swallowed by the Black Hole? If it could not be recovered, then 
the principle of unitarity of Quantum Mechanics would be violated. This problem is usu¬ 
ally known as the information paradox (see, for instance, [U] and references therein for a 
more accurate description of the problem). On the other hand, it is usually expected that 
quantum theory could avoid the existence of the singularity. Indeed, some paradigms and 
heuristic models of non-singular gravitational collapse and Black Hole evaporation inspired 
in different approaches to Quantum Gravity have appeared in the recent literature (see, for 
example, [Z][S]in][iD][ii][n][i3][ii]ii5]iis] and references therein). Even some specihc, more 
developed, models are now available (see, for instance, im for the collapse of shells, [IH] [IH] 
for homogeneous models or IZDl for scalar held collapse). 

Our aim in this work is to analyze the quantum corrected gravitational collapse of matter. 
We will do it by building and analyzing a family of models that tries to include some of the 
key physical phenomena involved in the process. In this way, the models will be able to be 
inhomogeneous and we will take into consideration that the generation of a horizon in the 
collapsing process requires the emission of Hawking radiation which will tend to evaporate 
the object. We will integrate the equations describing the collapse in order to analyze 
the complete stellar evolution. The hnal global outcome of gravitational collapse will be 
hypothesized accordingly, as well as the implications for the information paradox and for 
Ere walls. 

All the computations in this paper are carried out in Planck units where c = h = Gq = 1. 

The work has been divided as follows. In section we describe the general quantum 
corrected stellar interior. In section |3] we describe a stellar exterior suitable to describe 
Hawking radiation. The matching of the interior and the exterior solutions and a resolution 
scheme to obtain the full model can be found in section All the details for a particular 
model are worked out in section Finally, the results are discussed in the concluding section 
1^ where, based on our results, we conjecture a general picture for the global spacetime of a 
(non-singular) collapsing star that is able to generate a horizon. 
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2. Quantum corrected stellar interior 


In order to model the gravitational collapse of the quantum corrected matter we will assume 
that the whole spacetime will be split into two different regions V = V’*' U V“ with a common 
spherically symmetric time-like boundary E = dV^ fl dV~, corresponding to the surface of 
the star. We choose V’*' to represent the stellar exterior, while V~ will be its interior. In this 
section we will deal with the exact interior solution. 

Recently ini. it has been possible to provide the general solution for the quantum cor¬ 
rected collapse of non-interacting particles (dust) in a framework based on the Quantum 
Einstein Gravity approach to asymptotic safety. This class of ‘improved dust interiors’ is 
dehned by the metric 

dsi = -dr^ + ''}\ dr‘^ + R{t, Tfdn\ (2.1) 

1 -|- 2E{r) 

where r is the proper time of the particles composing the fluid, r is a parameter that labels 
every shell of the fluid, = dO"^ + seri^Odip'^, the prime indicates partial derivative with 
respect to r and R{T,r), the areal radius, is determined by 


if 

T 


GM{r) 

R 


E{r), 


( 2 . 2 ) 


with the overdot indicating partial derivative with respect to r, M(r) is the mass distribution, 
E{r) is the energy distribution and 


G = G{T,r) 


GoR^ 

R^+obGo{R + 'yGoM{r)) 


(2.3) 


is the running Newton constant. Here Go is Newton’s universal gravitational constant and 
Co and 7 are constants coming from the non-perturbative renormalization group theory and 
from an appropriate cutoff identihcation, respectively. The preferred value for these constant 
are 7 = 9/2 and Co = 167/307r. 

By choosing a function E{r) one chooses the total energy per unit mass of the particles 
in the fluid within a shell of radius r. Then, if E > 0 the system would be unbound, if E = 0 
the system would be marginally bound and if E < 0 the system would be doubly bound, 
meaning that there will not only be an upper bound to the particles, but also a non-zero 
(quantum) lower bound where the shell will bounce. 

The mass function [ 2 T[ | 22 l [23] of this solution is [H] 


M- = MG/Go 
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Figure 1: In order to illustrate the evolution of a given shell (i.e., a particular r) close 
to the center we have used the particular model that will he introduced in section 
Specihcally, we show the evolution of the areal radius R of the shell (solid line) 
as a function of its proper time. Likewise, we show the function R — 2GM along its 
evolution, what allows to see that starting from positive values of this function the 
exterior ASH is reached by this shell at r = 48.9519. Then, the shell traverses the 
region with trapped round spheres until it reaches the interior ASH at r = 51.9318. 


and the apparent 3-horizon (A3H) is the solution of the eqnation 

R - 2GM = 0. 

From all the ronnd spheres in the spacetime, only the ones satisfying R < 2GM are closed 
trapped snrfaces. 

The qnantnm corrected stellar dust behaves essentially as the classical dust far beyond 
the Planckian regime. In this way, starting from regular initial conditions, the model will 
collapse eventually generating an A3H. However, while in the classical case the later evolution 
of the shells will force them to fall into a singularity, the quantum corrected shells follow a 
different path: Once the shells cross the A3H they enter a region of trapped round spheres. 


While they go through this region the running G (2.3) decreases and, as a consequence. 


weakens the gravitational effects in such a way that, close to the Planckian regime, the shells 
cross an interior horizon. In this way, the shells enter an inner region of untrapped round 
spheres. No shell-focusing singularity forms [H]. We have illustrated this behaviour in hgure 
by showing the evolution of a particular shell in the particular model that will be used 
later in section [5l 

Let us now recall the physical content of the solution. Consider an observer with 4- 
velocity u = 8/dr and an orthonormal basis {u, n, such that ujq = R~^ 8/89, 

= (Rsm9)~^ 8/8(p and n is a space-like 4-vector. Supposing that it has been generated 
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by an effective matter fluid, the comoving observer with 4-velocity u will write the effective 
energy-momentum tensor for the anisotropic fluid as 


T =pu(8)u-hpn(8)n-Fj9_L {uJe ® uje + , (2.4) 

where p is the effective energy density, p is the effective normal pressure and p± is the effective 
tangential pressure. The held equations provide us with 

(GM)' _ MG,r , _ (MG,^)' 

^ ^ AnGoR SnGoRR' 


3. Stellar Exterior 

With regard to the stellar exterior region V^, we will describe it with a portion of a generalized 
Vaidya solution |25]. The spacetime metric for this solution can be written as 

dsl = - (l - du^ + 2dudR + Rdn\ (3.1) 

where dVR = dO'^ + sin^ Odp^ and Ad’*’ is the mass function of this solution. Depending on 
the specihc chosen mass function this metric has as particular solutions well-known classical 
solutions as the Schwarzschild solution (Ad'*' = constant), the Vaidya radiating solution 
(Ad"*" = Ad ■•■(«)) and others, but also quantum improved solutions such as the improved 
versions of the Schwarzschild solution (that implies a specihc dependence Ad"^ = Ad’''(.R)) 
[26] and the Vaidya solution (that implies a specihc dependence Ad"*" = Ad’''(M,.R)) [27][28]. 

The characteristics of this solution are very well suited for our model. On the one hand, 
we have seen that it is able to incorporate the required quantum corrections to the usual 
classical exterior solutions. On the other hand, this metric is able to describe the specihc 
evolving physical content of the stellar exterior: As the star begins its collapse in a phase 
previous to the formation of an ASH, we want the stellar exterior to be basically devoid of 
matter or radiation. This can be achieved if one considers for the exterior, at this initial stage, 
a portion of an improved Schwarzschild solution [2n]. Indeed, as intuitively expected, the 
model can be built by matching the improved interior and improved Schwarzschild solutions 
through a stellar surface (/matching hypersurface) S dehned by =constant [H]. However, 
once the collapse has reached the stage in which a horizon forms and the stellar surface is 
located inside it the existence of Hawking radiation in the stellar exterior implies that the 
exterior solution can no longer be an improved Schwarzschild solution. For example, there 
will be negative-energy fluxes (composed of radiation and, over time, of massive particles) 
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from the horizon falling towards S. Therefore, we should consider matching the interior 
improved solution with a particular portion of the generalized Vaidya solution possessing a 
more general dependence = M.'^iu^R). 

Note that the negative-energy fluxes could also traverse S and reach interior shells with 
a depth that depends on the translucency of the stellar interior. We will simplify the model 
by assuming that the star is perfectly opaque to the radiation so that the negative energy 
particles will just disintegrate the outermost shells of the collapsing star. In this way, the 
value of r on the matching hypersurface must decrease with time once the stellar surface has 
crossed the ASH. 


4. Matching surface with Hawking radiation 


In order to complete our model beyond the point the stellar surface reaches the ASH, let us 
now match the improved interior solution with the generalized Vaidya solution. V’''and V~ 
should satisfy Darmois matching conditions on E [29] [21], what implies that the solutions 
must be such that the hrst and second fundamental forms of S must coincide when computed 
from either or V~. 

We choose as a timelike parameter for the matching hypersurface the interior timelike 
coordinate r. The hypersurface will then be described, as seen from the interior, by r = r(r). 
The normalized 4-vector normal to E and the stellar surface’s 4-velocity can be written on 
E, from V“, as 


n = 


V = 


1 + 2E 


- 1/2 


d 1 + 2E d 




1 + 2E y / V 5 
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. d 

-L rp - 

dr dr 


1 + 2E -r^R'^^ 

As seen from V"*" {u = u{t), R = R{t)), they take the form, respectively. 


,R(i- 


2GoM- 

R~ 

2GoM- 


\ - 1/2 


— 2uR I ii -—h 


R 


\ - 1/2 


d 


4 1 - 


— 2uR ( ii -—h R-^ 




du 


d - d 


2GoM- 

R 


-R 


A 

M 


du 


dR 


(4.1) 

(4.2) 


(4.3) 

(4.4) 


Darmois matching conditions and, in particular, the requirement that the hrst funda¬ 
mental forms of E must coincide implies that the areal radii for the interior [R) and exterior 
regions {R) must agree on E [21]: 

R = R. (4.5) 
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Another consequence of the matching conditions are the Israel matching conditions 
[30] [23] that imply that the expression — n^) should take the same value when 

evaluated on E from the interior or the exterior. From the interior, one has 


T-0n^{v^ - n^) 


p{l + 2E — rR') + prR'irB! — 1) 

r 2/?'2 _ (1 + 2 E ) 


(4.6) 


while from the exterior 




{v^ - n^) = 


M~ 


yR 


47ri?2 


(4.7) 


If these two expressions must coincide on the matching hypersurface, then (using (4.5)) we 
obtain the evolution equation for the stellar surface. For example, in the case of an interior 
with E = 0 -we have 

s A-kR^p + 


r = 


R'{AttR‘^P-M- 




(4.8) 


4.1. Resolution scheme 


In order to model the collapse after E crosses the ASH we will proceed as follows. First, 
we will £x the stellar interior or, in other words, we will hnd a particular solution for 


(2.2). We have to choose particular energy {E{r)) and mass (M(r)) distributions for the 


interior together with an initial condition i?(ro,r). Then, ( 2.2[ ) will provide us a particular 
(numerical) solution i?(r, r). With the help of this solution one can ascertain the particular 


physical content of the interior (density and pressures) by computing (2.5). Second, we 
choose a particular exterior, i.e., a particular mass function Finally, the complete set 


of matching conditions provides us with (4.5), the evolution equation, the coincidence of 


mass functions at both sides of the matching hypersurface E 


M- =M- 


(4.9) 


and a differential equation for the evolution of the exterior time parameter u{t) (see [23] for 
details). The complete model can now be obtained from these conditions. In particular, one 
will now be able to write = E{u{t),R{t, r(r))) = E{t), so that the evolution equation 

may be numerically integrated. Incidentally, note that in the particular case of the evolution 


equation (4.8) a sufficient condition for the evaporation of the outermost shell of the star 


(r < 0) is F < 0 (since p > 0 and p < 0). 
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Figure 2: A plot of the areal coordinate R for shells around the center of the star, 
where one can follow the evolution of every shell (r=constant) according to its proper 
time T. As expected for E = 0 a the shells collapse towards the center. To plot 
this hgure we have arbitrarily chosen g = 1 and C = 10. 


5. A MODEL OF THE LAST STAGES OF STELLAR COLLAPSE 


We will choose for the interior a particular simple model in which E = 0, M{r) = gr and the 
initial condition is R{t = 0, r) = (r, where g and ( are constants. The numerical solution 
of (2.2) for a particular value of the constants provide us with the prohle for R{t,r) that 
can be seen in hgj^ The implementation of this particular solution in (2.5) provide us with 
the effective density and pressures that are shown in hgj^ As expected [H], note that the 
modulus of the density and pressures tend towards their maximum hnite Planckian valu^ 
given by 


pp ^ -Pp ^ -P±p ^ 


AnGlyu 


~ 0.02994. 


For our particular example in this section, we choose to work inside the ASH with a 
particular exterior solution that close to S takes the form|^ 


t>c 

^Note that the results are of the order of the Planck density and they are in Planck units. In SI units 
Pp ~ lO^^kg/rrd. 

^Similar and even nicer models can be implemented when the dependence of the exterior mass close to 
the total collapse is not linear in R but goes as i?" with n > 3. However, this does not affect significantly 
our main results. 
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Figure 3: In this figure we can check how the stellar interior evolves until reaching its 
maximum value for its density (pp ~ 0.02994j and its minimum value for the normal 
pressure (pp ~ —0.02994j. Note that p is always negative due to its quantum origin 

m- 


where Sq is a constant. (Incidentally, this solution has been recently used in String Theory 
EH in the context of radiating black holes). The solution has the advantage that it is 
relatively simple to work with for our specihc purposes. In particular, one can work the 
matching conditions in sequence. For instance, by using the previous numerical results for 
the stellar interior we can straightforwardly integrate the evolution equation for a given So 


. s SnR'^p — So 
R'{87rR^p + So)' 


(5.1) 


If the stellar surface must be a timelike hypersurface then, from (2.1) it has to satisfy 


-1 drs +1 

~W ^ ^ 

This requirement is clearly satished if one chooses a positive constant So- Moreover, this 
choice would also imply the evaporation of the star (r < 0). The result of the numerical 
integration of (5.1) is shown in hgj^ In the same hgure we also show the value of the function 
R — 2GM along the matching hypersurface in order to make explicit how the stellar surface 
traverses the region of closed trapped surfaces by crossing the inner apparent 3-horizon. The 


evolution of the function p with the internal time can be now obtained by using (4.9). This 
provide us with 


p[t) = /i(r, r(r)) = M(r(r))G(r, r(r)) + i?(r, r(r))^o/2. 

This is a decreasing function with time that reaches zero when the stellar surface completely 
collapses (i? —0, what implies G —)■ 0). 
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Figure 4: The evolution of the areal radius of the stellar surface Rt,{t) (solid line) has 
been obtained once it crosses the exterior ASH. The value of the function R — 2GM 
along S is also shown (dashed line) as a reference. In this particular case, the stellar 
surface leaves the region of trapped round spheres when it crosses the interior horizon 
at T = 50.5224. The total evaporation event is at t = 52.3912, when S reaches R = 0. 
The specific initial condition on the ASH is the one in the caption of fig\^ and we 
have arbitrarily chosen the particular value Sq = 0.01. 

6 . Conclusions 

In this paper we have addressed the problem of the quantum corrected gravitational collapse 
of a star. The modeled stars had an interior composed of non-interacting particles where 
quantum gravity effects had been taken into consideration (Secj^. These quantum effects 
make the interior free of shell-focusing singularities. The interior has been matched to an 
exterior that is suitable to describe quantum corrections outside the star as well as the 
Hawking radiation that will be produced once an ASH is generated (Secj^. We have taken 
into account that the negative energy flux of Hawking radiation that will then reach the 
(opaque) star should evaporate its outermost shells. Therefore, the stellar surface S that is 
hrst described by =constant satishes, once Hawking radiation reaches it, rs < 0 (Secj^. 
In this way, a total evaporation (r^ = 0) due to Hawking radiation seems feasible. In 
principle, the process would satisfy energy conservation since the initial content of the star 
that is consumed by the negative energy flux of Hawking radiation is compensated by an 
equivalent outward flux of positive energy photons and massive particles. 

In order to illustrate these ideas we have considered a particular conceptual model (Secj^ 
and we have analyzed it close to the total evaporation event. The interior has been chosen 
so that shell-crossing singularities are avoided in the region under consideration, what is 
the usual approach since possible shell-crossing singularities are just mathematical artifacts 
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caused by dealing with noninteracting dust shells. Regrettably, in order to the shells bouncing 
at the center not to cross between them, this forces us to choose a specihc model that quickly 
collapse^ This is not totally satisfactory since one would expect the shells crossing the 
inner horizon to coexist (an interact) for a longer time while trapped by the inner horizon 
(0 < i? < Ria 3 h)- This trap should eventually provide the collapsed star with enough time 
to evaporate as all the generated negative flux of Hawking radiation reaches it. Therefore, 
a more complete model would require the description of interacting shells. However, this 
possibility is not available at the present state of the art and it should be left for future 
works. 

The global picture of gravitational collapse that one could extract from our general and 
particular models, as well as our arguments above, can be seen in £gure[^ It is important to 
emphasize that the global structure of the spacetime is devoid of an event horizon. Therefore, 
according to the standard dehnition there is no black hole. However, from the point of view 
of some external observers the object that forms is not distinguishable from a standard 
evaporating black hole during the evaporation time (about 10^*^ years for a solar mass black 
hole). The collapse is completely regular since no singularity appears. The information 
loss problem does not have any meaning since there is not a singularity destroying the 
information, the apparent 3-horizon (that replaces the event horizon) can be traversed by 
both ingoing and outgoing particles and the spacetime global structure allows the information 
about the collapsed matter to scape towards null inhnity at late times. Note that there is 
no need for a firewall in the vicinity of the ASH in order to preserve unitarity. In this way, 
we have obtained a picture for gravitational collapse that seems reasonable and satisfactory. 

Only a warning must be considered: It has been argued |26] |T3] that the quantum correc¬ 
tions used in this work are accurate for R of the order of Planck’s length. Therefore, despite 
we have depicted the full evolution of the star until its complete evaporation, we cannot 
guarantee that the last stages of the collapse will be exactly as depicted. The Planckian 
region where only a full theory of Quantum Gravity could make a totally reliable prediction 
of the events has been emphasized in hgj^ 

An interesting feature of our model for the global spacetime is that there are external 
observers in the causal future of the Quantum Gravity domain (£g|^. Not only will these 
observers be able to recover the apparently lost correlations of matter and radiation that 
had fallen into the ASH, but they will also be able to collect experimental data from the 
Planckian region. In this way, in principle, it could be possible to test different theories of 
Quantum Gravity by observing evaporating Black Holes. 

®I.e., shell-crossing is clearly unavoidable in this model if we had allowed the innermost shells to bounce 
at the center. 


11 



R=0 



Figure 5: The Penrose diagram of gravitational collapse that one could infer from 
the computations in this work. The stellar interior is limited by the stellar surface 
(matching hypersurface) S that collapses towards R = 0. In the first stages the 
exterior in basically vacuum and =constant. An apparent 3-horizon (ASH) even¬ 
tually forms and, in this quantum corrected case, it is closed and contains a region 
of spherically symmetric closed surfaces with R > 0 (stripped region) so that one 
can distinguish between an interior ASH and an exterior ASH. The formation of the 
ASH implies the emission of ingoing and outgoing Hawking radiation in the form of 
photons and massive particles. We have just depicted the outgoing photons by using 
wavy arrows. The rest of the generated radiation and the remains of the evaporating 
star are symbolized by a fading grey in the exterior region. The dashed line around 
the complete evaporation event represents the boundary of the planckian sized region 
which should be better described by a full theory of Quantum Gravity. 
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